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ABSTRACT 



We calculate the radial surface brightness profile of the image of a Gamma- 



Ray-Burst (GRB) afterglow. The afterglow spectrum consists of several power-law 
segments separated by breaks. The image profile changes considerably across each of 
the spectral breaks. It also depends on the density profile of the ambient medium into 
which the GRB fireball propagates. Gravitational microlensing by an intervening star 
■ can resolve the afterglow image. We calculate the predicted magnification history of 

GRB afterglows as a function of observed frequency and ambient medium properties. 
{Sj , We find that intensive monitoring of a microlensed afterglow lightcurve can be used to 

reconstruct the parameters of the fireball and its environment, and provide constraints 
£T) ' on particle acceleration and magnetic field amplification in relativistic blast waves. 
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Introduction 



45 

Or 
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The fireball of a Gamma-Ray Burst (GRB) afterglow is predicted to appear on the sky 
as a ring (in the optical band) or a disk (at low radio frequencies) that expands laterally at a 
superluminal speed, ~ Tc, where r ^> 1 is the Lorentz factor of the relativistic blast wave which 
emits the afterglow radiation (Waxman 1997; Sari 1998; Panaitescu & Meszaros 1998; Granot, 
Piran, k. Sari 1999a,b). Days after the GRB trigger, the physical radius of the afterglow image 
(~ the fireball radius over T) translates to an angular size of order a micro-arcsecond (//as) at a 
cosmological distance. Coincidentally, this image size is comparable to the Einstein angle of a 
solar mass lens at a cosmological distance, 

where Mi ens is the lens mass, D = D 0S D \/ D\ s is the ratio of the angular-diameter distances 
between the observer and the source, the observer and the lens, and the lens and the source 
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(Schneider, Ehlers, & Falco 1992). Loeb &i Perna (1998) argued that microlensing by stars could 
therefore be used to resolve the photospheres of GRB afterglows. 

Recently, Garnavich, Loeb, & Stanek (2000) have reported the possible detection of a 
microlensing magnification feature in the optical- infrared light curve of GRB 000301C. The 
achromatic transient feature is well fitted by a microlensing event of a 0.5M Q lens separated by an 
Einstein angle from the source center, and resembles the prediction of Loeb & Perna (1998) for a 
ring-like source image with a narrow fractional width (~ 10%). Alternative interpretations relate 
the transient achromatic brightening to a higher density clump into which the fireball propagates 
(Berger et al. 2000), or to a refreshment of the decelerating shock either by a shell which catches 
up with it from behind or by continuous energy injection from the source (e.g., Zhang & Meszaros 
2000). In order to regard these alternatives as less plausible, it is essential to identify the unique 
temporal and spectral characteristics of a microlensing event. 

Gravitational lensing of a point source is achromatic due to the equivalence principle. 
However, differential magnification by an achromatic lens could still produce a chromatic 
magnification signal if the source is extended and its image looks different at different photon 
frequencies (see, e.g. Loeb & Sasselov 1995). Should a microlensing event of GRB afterglows 



be achromatic? Since the image profile of an afterglow does depend on frequency (Sari 1998 



Panaitescu fc Meszaros 1998| ; Granot, Piran, & Sari 1999a, b), its magnification lightcurves must 
exhibit a predictable level of chromaticity. In §2 we will demonstrate that the surface brightness 
profile of an afterglow image changes considerably with frequency across all spectral breaks in its 
broken power-law spectrum. We will therefore show in §3 that the magnification history of an 
intervening microlens would depend on frequency across spectral breaks. However, within each 
power-law segment of the afterglow spectrum, the magnification would remain achromatic, in 
rough consistency with the sparse optical-infrared observations of GRB 000301C. This spectral 
behavior offers a fingerprint that can be used to identify a microlensing event and distinguish 
it from alternative interpretations. It can also be used to constrain the relativistic dynamics of 
the fireball and the properties of its gaseous environment, as well as the length scale required for 
particle acceleration and magnetic field amplification in relativistic shocks. 



2. Frequency Dependence of Afterglow Images 

The surface brightness profile (SBP) of an afterglow depends on the underlying hydrodynamics 
of the fireball, as well as on the energy distribution of shock-accelerated electrons. In this Letter 
we adopt the model used by Granot & Sari (2001, hereafter GS01) to describe afterglows. Here 
we provide a brief outline of this model, and refer the reader to Granot et al. (1999a,b) and GS01 
for more details. 

The hydrodynamics is described by the Blandford-McKee (1976) self-similar spherical solution, 
with a power-law external density profile, p cx R~ k , for either k = or k = 2, corresponding to a 
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uniform interstellar medium (ISM) or a stellar wind environment, respectively. The number per 
unit energy of accelerated electrons is assumed to be a power law, dN e /d^ oc 7~ p (for 7 > ~f m ), 
just behind the shock. Thereafter, the initial distribution evolves due to radiative and adiabatic 
losses. The emission mechanism is assumed to be synchrotron radiation, and the emissivity is 
integrated over the entire volume behind the shock front. 

The afterglow image is limited to a circle on the sky, whose size grows as 
is the observed time. The assumption of a spherical flow may also serve as an adequate description 
of a jetted flow, at sufficiently early times before the jet break time, tj e t, when the Lorentz factor 
of the flow drops below the inverse of the jet opening angle (Rhoads 1997). In the case of a jet, 
the image at t ^ tj e t is expected to be different than in the spherical case, and will no longer be 
circular for observers who are situated away from the jet axis. Deviations from sphericity might 
also result from hydrodynamic or plasma instabilities. For simplicity we focus here on a spherical 
fireball; the more complex cases mentioned above are left for a future study. 

The spectrum of GRB afterglows consists of several power-law segments (PLSs) where the 
flux density F v oc v@, which join at certain break frequencies ( |Sari, Piran Narayan 199S , 



hereafter SPN98; Pranot, Piran E Sari 2000| , hereafter GPSOO; GS01). The cooling time of a 
typical electron, t coo i, is initially smaller than the dynamical time, t<fyn> resulting in "fast cooling" 
(SPN98), while at latter times t coo i becomes larger than td yn and there is a transition to "slow 
cooling". The fast cooling stage is typically expected to last for the first hour (day) for an ISM 



(stellar wind) environment (SPN98; Chevalier &: Li 2000 ). Since the angular size of the source is 



very small at such early times, the source will remain <C Be and hence unresolved during the fast 
cooling regime of a typical microlensing event (or else be resolved only in a negligible fraction of 
all events for which the lens nearly coincides with the source center) . We will therefore restrict our 
attention to the slow cooling regime. The slow cooling spectrum may assume one of three different 
shapes, depending on the ordering of the self-absorption frequency, u a , with respect to u m < u c , 
where v m is the typical synchrotron frequency and v c is the cooling frequency (GS01). The most 
common ordering of the break frequencies is the one shown in Figure [I] with v a < u m < v c . For 
clarity, we show a broken power law spectrum, even though in the actual spectrum obtained from 
our model the PLSs join smoothly, and the break frequencies are defined where the asymptotic 
power laws meet (e.g. Granot et al. 1999b). Altogether, there are five different PLSs that appear 
in the different slow cooling spectra, which correspond to /3 = 5/2, 2, 1/3, (1 — p)/2, —p/2. Due to 
the self-similar hydrodynamics, the surface brightness, normalized by its average value over the 
image, is independent of time within a given PLSg. For (3 = {l-p)/2 and (3 = -p/2, the value of 
(3 depends on p. The values of (3 in these two PLSs may therefore accidentally coincide with those 
of other PLSs. However, the SBP will still be different, since the underlying physics is different. 



3 This statement remains valid even when taking into account the fast cooling spectra, which include two additional 
PLSs corresponding to = 11/8, -1/2 (e.g. SPN98; GPSOO), with only one exception: (3 = 1/3, where the SBP is 
different between slow cooling and fast cooling, since the underlying physics is different. 
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Although for these two PLSs the SBP depends on the value of p, we find that at any given time, 
the differences in the microlens magnification induced by the differences in the SBP, never exceed 
a few percent for 2 < p < 3 (and b ~ 1, see §[|). Hence, we use a single value of p = 2.5 in the 
following. 

Figure || shows the surface brightness, normalized by its average value, as a function of the 
normalized radius inside the image, r, for k = 0, 2 and (3 = 5/2, 2, 1/3, (1 — p)/2, —p/2. The SBPs 
for k = and [3 = 2, 1/3, (1 — p)/2 are taken from Granot et al. (1999a, b), while for the other 
values of k and /3, the SBPs are calculated in a similar manner, using the formalism of GS01. 
For both values of k, the image tends to be more uniform as (3 increases and more ring-like (i.e. 
dim near the center and bright near the outer edge) at low values of (3. For a given PLS (i.e. a 
given (3 in Fig. the image tends to be more uniform for a stellar wind environment (k = 2) as 
compared to an ISM environment (k = 0). For (3 = —p/2 the surface brightness diverges at the 
outer edge of the image (r = 1), since the emission in this case is both optically thin and originates 
from a very thin layer behind the shock, leading to a relativistic version of the limb brightening 
effect (e.g. Sari 1998| ). However, the physical size of this thin layer cannot be smaller than that 



of the transition layer just behind the shock front, where the electrons are accelerated and the 
magnetic field is amplified. As will be shown in §||, the divergence of the SBP at r = 1 leads to 
a very sharp spike in the magnification history. Introducing a finite width to the transition layer 
would eliminate the divergence of the SBP. However, a thin transition layer would still result in 
a sharp peak of the SBP near r = 1, and a sharp peak in the magnification history, whereas a 
thick transition layer will result in a lower and smoother peak in both. Therefore, the sharpness of 
the magnification spike may set constraints on the physical processes of particle acceleration and 
magnetic field amplification in relativistic blast waves. 



3. Microlensing Histories 

Next we calculate the magnification histories of GRB afterglows due to microlensing by an 
intervening star at a cosmological distance. Since the effects of microlensing depend on the angular 
structure of the image in units of the Einstein angle Be [Eq. (|lj)], we normalize all angular scales 
by #e- The angular radius of the afterglow image, in units of 0Ej is given by 

Rsit) = R,tt: )m - k \ (2) 

where tdays is the observed time in days and Rq is the angular size of the image in units of #e after 
one day. The angular separation between the center of the image and the lens, in units of #Ej is 



4 The divergence scales as (1 — r) _1//2 , since the intensity from an optically-thin shell is proportional to its geometric 
length along the direction to the observer. This limb brightening effect does not occur for the other optically-thin 
PLSs we consider, where the emission originates from the entire three-dimensional volume behind the shock. Only 
the fast-cooling electrons produce the radiation from a very thin layer just behind the shock front. Note that the 
diverging cusp makes a vanishing contribution to the average surface brightness of the afterglow image. 
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denoted by b. The magnification of a uniform ring of fractional width W (occupying the interval 



1 — W < r < 1) is given by (Loeb & Perna 199 



KRs, W, b) _____ ( 3 ) 



where *&(R S , b) is the magnification of a uniform disk of radius R s (Schneider et al. 1992), 

r b+R s D2 , 9 A2 , r>2 _ t>2 ' 



*(Rs, b) 



ttR 2 



R ° R 2 + 2 b 2 + R 2 - R 2 

e - + H{K s -b)-{R s -b)^{R s -b)-+4 

(4) 

and H(x) is the Heaviside step function. In order to obtain the magnification of a radially 
symmetric image, one can either divide it into a large number of thin rings, or equivalently use the 
resulting integral formula for the frequency-dependent magnification (Witt Sz Mao 1994), 



Hu(Rs, b) = [ 2rdr 
Jo 



*(rR s , b) + ~^(rR s , b) 



t(r) (5) 



where l v = I V /{I V ) is the specific intensity normalized by its average value, = Jq 1 2rdrl u (r), as 
shown in Figure |2|. Note that for a uniform disk (I v {t) = 1), equation (|5|) gives \x v = ^f(R s , b), as 
it should. 

The magnification history for b = 1, k = 0, 2 and (3 = 5/2, 2, 1/3, (1 — p) /2, —p/2, is shown in 
Figure ||. The peak in the magnification history is generally higher and sharper, and occurs at 
a slightly earlier time for a smaller spectral slope, (3. This follows from the general trend of the 
SBP, which tends to be more sharply peaked near the outer edge of the image for lower values 
of (5. At high values of (3, especially in the optically-thick regime (f3 = 5/2,2), the magnification 
peak is relatively shallow and broad. For an ISM environment (k = 0), the peak magnification is 
somewhat higher and the differences in the magnification histories for different (5 values are more 
pronounced than for a stellar wind environment (k = 2). These differences may allow to uncover 
the density profile of the circumburst medium from intense monitoring of the microlensing 'bump' 
in an afterglow light curve. Unfortunately, the quality of the observational data for GRB 000301C 
is not sufficient for this purpose. 

We now consider the effects of microlensing on the global spectrum. For clarity, this is 
illustrated in Figure |l], on the spectrum therein. The spectrum consists of several PLSs: (3\ = 2, 
/?2 = 1/3, 03 = (1 — p)/2 and (3% = —p/2, with different magnification histories, fii(t) (i = 1, ...,4). 
This causes each PLS to shift by a different time-dependent factor, as illustrated by the thin solid 
line in Figure [l] for t^ays = 1- The different magnification in each PLS generates a shift in the 
location of the break frequencies compared to their unlensed values, as illustrated in the insert at 
the top right corner of Figure |l[ The value of the z-th break frequency (numbered from low to high 
frequencies, as are the PLSs) will shift by a factor of [/x i+ i(t)//Xi(t)] 1/(ft_ft+l) . The largest shift 
occurs at v c since the change in the spectral slope across the break is rather small {(3% — [3^ = 1/2) 
and the ratio [ii{t) / '^(t) reaches relatively high values (see the upper panel of Figure H). 
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4. Discussion 

We have shown that the frequency dependence of the afterglow image profile (Fig. 2) 
introduces a specific level of chromaticity to the magnification history of GRB afterglows (Fig. 
3). This unique fingerprint can be used to distinguish between microlensing and alternative 
hydrodynamic interpretations for the transient brightening in an afterglow light curve. 

For example, if the GRB blast wave runs into a clump of gas along the line-of-sight (Berger et 
al. 2000), then brightening will occur first in the radio and later in the optical, in opposite order 
to the microlensing case. The chromaticity in the clump case is sensitive to the exact geometry of 
the clump and the sharpness of its boundary (which could send a reverse shock into the fireball). 
In another example, the brightening event may result from refreshment of the decelerating shock 
(Zhang & Meszaros 2000). In this case, the break frequencies and peak flux of the spectrum in 
Figure |l|, namely v a , u m , v c and F Um , will change by a factor of f 1 ^, f 1 ^ 2 , / _1//2 and / (or by 
/~ 2/ ' 5 , / 1//2 5 / 1//2 and J 1 / 2 ), respectively, for an ISM with k = (or stellar wind with k = 2) 
environment, where / is the ratio between the total energies in the forward shock after and before 
the refreshment. This behavior clearly differs from the microlensing prediction (see insert in Fig. 
H). By itself, the magnification history of microlensing is already unique in that it is initially 
characterized by a constant positive offset as long as R s (t) <C 1 (Loeb & Perna 1998). Such an 
initial positive offset does not appear naturally in either of these alternative interpretations. 

Microlensing events are rare but precious. Only one out of roughly a hundred afterglows is 
expected to be strongly microlensed (Press h Gunn 1973; Blaes & Webster 1992; Koopmans & 
Wambsganss 2000), although all afterglows are expected to be magnified at some weak level (Mao 
&c Loeb 2000). Detailed monitoring of a few strong microlensing events among the hundreds of 
afterglows detected per year by the forthcoming Swift satellite^, could be used to constrain the 
environment and the dynamics of relativistic GRB fireballs, as well as their magnetic structure 
and particle acceleration process (Fig. 3). Unfortunately, the data on the existing microlensing 
candidate GRB 000301C is not of sufficient quality to provide firm constraints on these parameters, 
although the i?-band data seems to suggest (3 = —p/2 and gives a somewhat better fit for k = 
(a uniform density environment) than for other parameter choices. This can be seen by crudely 
approximating the SBP as a ring. We find that for k = 0, a fraction of 2/3 (3/4) of the emission 
comes from a ring of fractional width W = 0.156 (0.190) for = (1 - p)/2 and W = 0.075 (0.104) 
for /3 = -p/2; while for k = 2, W = 0.318 (0.370) for (3 = (1 -p)/2 and W = 0.128 (0.176) for 
(5 = —p/2. The k = 0,f3 = —p/2 values get closest to the empirical best-fit values of W = 0.07 
using the R-band data (Garnavich et al. 2000). The statistical accuracy by which the unknown 
brightness profile of a future microlensed afterglow can be reconstructed through an intensive 
monitoring campaign, was analyzed in detail recently by Gaudi & Loeb (2001). 



5 Planned for launch in 2003; see 



http:/ /swift. sonoma.edu/ 
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Fig. 1. — Typical broken power-law spectrum of a GRB afterglow at a redshift z = 1 (boldface solid 
line). The observed flux density, F u , as a function of frequency v, is shown at an observed time 
tdays = 1 f° r an explosion with a total energy output of 10 52 erg in a uniform ISM (k = 0) with a 
hydrogen density of 1 cm -3 , and post-shock energy fractions in accelerated electrons and magnetic 
field of e e = 0.1 and e# = 0.03, respectively (using the scalings and notations of GS01). The thin 
solid line shows the same spectrum, when it is microlensed by an intervening star with 6=1 and 
Rq = 1 (see §^). The insert shows the excess evolution of the break frequencies fb = v a , v m and v c 
(normalized by their unlensed values) due to microlensing. 
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Fig. 2. — The surface brightness, normalized by its average value, as a function of the normalized 
radius, r, from the center of the image (where r = at the center and r = 1 at the outer edge). The 
image profile changes considerably between different power-law segments of the afterglow spectrum, 
F u oc v 13 . There is also a strong dependence on the density profile of the external medium, p oc R~ k . 
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Fig. 3. — The magnification of the afterglow, fi, as a function of normalized time, for different 
power-law segments, F v oc v^, and different external density profiles (p oc Rr k for k = 0,2). The 
lens is assumed to be separated by one Einstein angle from the center of the source (b = 1). 



